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Determining the optimal arrangement of superconducting layers to withstand large amplitude 
AC magnetic helds is important for certain applications such as superconducting radiofrequency 
cavities. In this paper, we evaluate the shielding potential of the superconducting film/insulating 
hlm/superconductor (SIS’) structure, a configuration that could provide benefits in screening large 
AC magnetic helds. After establishing that for high frequency magnetic helds, hux penetration must 
be avoided, the superheating held of the structure is calculated in the London limit both numerically 
and, for thin hlms, analytically. For intermediate him thicknesses and realistic material parameters 
we also solve numerically the Ginzburg-Landau equations. It is shown that a small enhancement of 
the superheating held is possible, on the order of a few percent, for the SIS’ structure relative to a 
bulk superconductor of the him material, if the materials and thicknesses are chosen appropriately. 

PACS numbers: 74.78.Fk, 74.25.0p, 64.60.Ht, 77.55.-g 


I. INTRODUCTION 

Can one engineer a better superconducting magnetic 
shield? How can one optimally arrange materials to 
maintain complete flux exclusion from a region, and what 
is the maximum external field that can be screened? It 
has long been known that superconducting films of width 
d smaller than the London magnetic penetration depth 
A can remain superconducting at much higher magnetic 
fields than bulk samples,so it has been proposed that 
films could be used to shield bulk superconductors.^ In 
this paper, we investigate the shielding properties of the 
film/insulator/bulk (SIS’) structure and compare to the 
single superconducting slab. The focus here is on AC 
rather than DC shielding; the latter has already been 
studied extensively.l^lt^ 

Superconducting radio-frequency (SRF) cavities are an 
example of an application in which shielding of large- 
amplitude high frequency magnetic fields is required. 
This technology underlies particle accelerators used in 
high-energy physics, nuclear physics, neutron sources, 
and X-ray light sources. The large AC accelerating elec¬ 
tric field of these cavities induces a correspondingly large 
magnetic field. If the magnetic field exceeds the flux pen¬ 
etration field of the material, it causes a quench in the 
cavity. If SIS’ structures could enhance the flux pene¬ 
tration field relative to that of a bulk superconductor, it 
could allow these cavities to achieve higher accelerating 
fields.^ In this paper, we examine the superheating fields 
Bsh of these structures, where flux penetration would oc¬ 
cur in defect-free superconductors. In fact, part of the 
motivation for this work is that there has been significant 
confusion in the SRF community regarding the maximum 
fields that SIS’ structures can screen; we hope that this 
study will clarify the screening mechanism and its limita¬ 
tions. Our calculations show modest shielding gains for 
SIS’ heterolaminates compared to bulk superconductors. 


The SIS’ structure may provide benefits in other ways for 
realistic materials with surface defects,!^ but considering 
those benefits is beyond the scope of the present work. 

The paper is organized as follows: we start our analysis 
by arguing in Sec. |TT] that for an SIS’ structure, a signif¬ 
icant enhancement of the flux penetration field could be 
achieved only if a significant gradient in the phase of the 
order parameter V/ can be established across the film 
shielding the bulk. Since this would result in a level of 
dissipation that is likely unmanageable, we restrict our 
analysis to fields below Bgh, where both film and bulk 
superconductor are in a (meta)stable regime, with no 
phase gradient across the film. In Sec. [ml numerical cal¬ 
culations are performed in the London limit. The thin 
film regime is examined in Sec. |IV| with an analytical 
Ginzburg-Landau approach. In Sec. |Vj the results are 
extended to Aims of intermediate thicknesses via a full 
numerical Ginzburg-Landau analysis. We summarize our 
work in Sec. IVII 


II. FLUX EXCLUSION 

The fundamental link between superconducting order 
and magnetism is the fact that the free energy and prop¬ 
erties of the system are governed not by the gradients Wtjj 
of the superconducting order, but by a ‘covariant’ deriva¬ 
tive Dij) = (V — e*iA/h)ip, where e* = 2e is the Cooper 
pair charge and A is the magnetic vector potential. If we 
write the complex superconducting order parameter in 
terms of two real fields as t/j = \ijj\ exp(i())), the covariant 
derivative becomes 

Dij} = [V|/| -I- i|/|(V(() — e*A/h)] exp(i/) (1) 
= [V|'0| -b \\%|^\{m*^f,/h)] exp(i/), (2) 
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where the gauge-invariant combination 

— e* A/h) = Vs (3) 

is called the supercurrent velocity. Magnetic fields cause 
“stress” in superconductors indirectly through A, which 
induces screening supercurrents. Due to these supercur¬ 
rents, a weak magnetic field exponentially decays inside 
a superconductor over the penetration depth A. As a 
crude approximation, the superconductor can support a 
certain maximum stress, characterized by a maximum 
superfiuid velocity The superconductor can only 

screen A values larger than je* if it passes vor¬ 

tex lines through its boundary. For example, if a vortex 
line is passed through a hollow superconducting cylinder 
in a parallel external field, this will bring flux inside the 
cylinder and (j) will wind by a factor of 27r, lowering the 
stress in the superconductor. 



X [nm] 


through the film, the vortices would experience strongly 
dissipative drag,^ generating levels of heating that are 
likely unmanageable for most applications 

As we are focusing on high frequency applications, we 
impose the restriction that flux must never pass into the 
superconducting regions. With this restriction, the SIS’ 
structure would offer an advantage over a single thick 
superconducting slab if it could withstand higher mag¬ 
netic fields without flux penetration. The flux-free state 
is only intrinsically stable below i?ci, the lower critical 
field. However, there is good evidence that real mat erials 
can withstand high frequency fields well above 
In this metastable regime, an energy barrier prevents flux 
from penetrating, a barrier that is reduced to zero at Bgh 
for a defect-free material (thermal fluctuations at cryo¬ 
genic temperatures are much smaller than the condensa¬ 
tion energy, so they cannot create excitations above the 
barrier). Bsh is the ultimate AC magnetic limit; this is 
especially important for SIS’ films, as they are always in 
the metastable state.l^ We will use B^^^ to denote the 
maximum metastable field of a SIS’ structure to distin¬ 
guish it from the superheating field of the bulk material, 
Bsh,b, and the bulk superheating field of the film mate¬ 
rial, Bshj (i-e. the value it would have if it were not 
a thin film). In the next three sections we present and 
compare three approaches to evaluate . 

III. SUPERHEATING FIELD IN THE LONDON 
LIMIT 


FIG. 1: Example of a SIS’ structure. The amplitudes of the 
magnetic field, the vector potential, and the Gibbs free energy 
are plotted as a function of distance into the structure. 

Now let us consider a single thin superconducting film 
separated from a bulk superconductor by a thin insu¬ 
lator, shown in Figure In a ‘thin’ superconductor of 
thickness d X the critical fields are enhanced; for exam¬ 
ple, for the parallel thermodynamic critical field we have 
iJcii = The Meissner state requires A —0 

deep in the bulk, and A is continuous across the insu¬ 
lating gap. Therefore, the vector potential at the film 
surface is tied to that of the bulk superconductor sur¬ 
face; however, the insulating gap offers the opportunity 
to decouple the phase gradient across the film from that 
in the bulk. If many vortex lines pass through the film, 
the superconducting film could be relatively unstressed, 
supplementing the native superheating field of the film 
material. 

In DC, it should be possible to screen a bulk from very 
large fields using a compound film with many layers of al¬ 
ternating thin superconducting and insulating films with 
magnetic flux trapped between each of them. However, in 
AC applications, filling the insulators with magnetic flux 
demands the transfer of V())/7r fluxoids per unit length 
across the screening film in each cycle. As they pass 


To make a rough estimate of the superheating field of 
the SIS’ structure, we consider the Gibbs free energy Q in 
the London limit; that is, we assume that both film and 
bulk superconductors are strongly type H materials, with 
penetration depths much longer than coherence lengths. 
We denote by A/ the film’s material penetration depth 
and by its coherence length. The thickness d of the 
film is assumed to be much larger than in particular, 
for the vortex core t o be accommodated in the film one 
needs d > The film is separated from a bulk 

superconductor with with penetration depth A{, by an 
insulating film of thickness 5. The superconducting film 
is screening the bulk from a parallel magnetic field with 
amplitude Bq. The screened field between the film and 
the bulk has amplitude B^. In our geometry, the x-axis 
is perpendicular to the film, pointing into it, with origin 
at the interface with the exterior. The 2 :-axis is aligned 
with the magnetic field. 

The Gibbs free energy of a vortex in a superconduc¬ 
tor can be determined from the value of two magnetic 
fields evaluated at the vortex location tq: the Meissner- 
screened external field Bm and the field generated by the 
vortex in the film HylSl 

g=^{Bv{ro)/2 + BMiro)) , (4) 

Mo 

where (po is the flux quantum and po the magnetic con- 
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Material 

A [nm] 

e [nm] 

Bel [T] 

B. [T] 

Bsh [T] 

Nb 

40 

27 

0.13 

0.21 

0.25 

NbgSn 

111 

4.2 

0.042 

0.50 

0.42 

NbN 

375 

2.9 

0.006 

0.21 

0.17 

MgB2 

185 

4.9 

0.017 

0.26 

0.21 


TABLE I: Materials parameters of niobium and three 
promising alternative SRF materials. The penetration depth 
A is calculated using Eq. 3.131 in Ref. 1151 The correlation 
length ^ is calculated using the equations in Ref. 1161 For 
Nb a RRR of 100 was assumed. For MgB 2 , A and ^ are not 
calculated, as the experimental values are given in the ref¬ 
erence. For calculations, = 0o/(2\/27r^A) is usedRei 
for Nb found from power law fit to numerically computed 
data from Ref. IT] and and for strongly type II materi¬ 
als is found from Eq. 5.18 in Ref. HU Bsh is calculated using 
Bsh — Be from Ref. 1191 Nb data from Ref. 1201 

NbsSn data from Ref. 1211 NbN data from Ref. 1221 and MgB 2 
data from Ref. 1231 Note that the two gap nature of MgB 2 
may require more careful analysis than is performed here. 



FIG. 2: Maximum field below B^h of both the film and the 
bulk as a function of film thickness for various film materials in 
a SIS structure with Nb. The effect of varying the insulator 
thickness 5 is shown for the NbaSn film, as is the effect of 
splitting the film thickness d over 5 equally thick multilayers 
with thin separating insulators. All calculations done in the 
London limit. 


stant. The field Bm can be found by minimizing the free 
energy in the structure when no vortex is present; we re¬ 
mind that in the London limit the Meissner field in the 
bulk superconductor decays exponentially and hence it 
equals . This procedure gives: 


Bm = 


1 x—dl2 • -I x — dl2 

Bn + B, cosh Bn - B, sinh 


cosh ^ 2 sinh ^ ’ 

2\f 2\f 


where Bi is given by 
Bi = Bq 


S + \b . d d 

—r sinh h cosh — 


( 5 ) 


( 6 ) 


Explicit formulas for By are available for thin {d <C 
A/) and thick {d^ A/) films.^i^To study the full range of 
thicknesses, we use the more general expression of R ef. l24l 
(this expression assumes tq = (xo,0)): 


The structures plotted in Figure can be divided into 
two types: homolaminates, in which the film is the same 
material as the bulk, and heterolaminates, in which they 
are different. Calculations show that for a homolami¬ 
nate like NbaSn/insulator/NbaSn, the film is the weak 
point: it always reaches its Bgh before the bulk, and the 
thinner the film, the lower its B^h- Homolaminates with 
films that are so thick that they behave like a bulk super¬ 
conductor have the highest . To better understand 

this, consider the magnetic forces on a vortex [which can 
be derived from Eq. @]. The boundary condition im¬ 
posed by By can be satisfied by an image antivortex 
outside of the boundary, which creates a force that pulls 
the vortex out of the film.l^ As the film thickness is re¬ 
duced, Bm remains approximately unchanged, but the 
image antivortex on the insulator side of the film used to 
satisfy By has a stronger effect, as shown in Figure 
This lowers the barrier to penetration. 


f dk ifct/ sin(7rna:/ d) sm{7rnxo/d) 

^ X^d^ J 2^® fc2 + (^n/d)2 + l/A2 

— oo 

Equations @-0 give the fields in the structure, and 
Eq. Q gives the Gibbs free energy as shown in Fig. 
The barrier to flux penetration is due to the positive slope 
of Q inside the superconducting regions near the inter¬ 
faces. We can find Bf^^ by finding the field at which the 
barrier is reduced to zero in any of the superconductor^^. 
In Fig. i BfA is plotted as a function of superconduct¬ 
ing film thickness for various SIS’ structures. Various 
insulator thicknesses are considered, including the thin 
layer limit, for illustrative purposes as it gives the highest 
fields. The materials analyzed are those that are promis¬ 
ing for SRF cavities, with properties given in Table |T] 


Meissner-screened external 



Exterior Thin Film SC Insulator 



Image 
antivortex 
V (attracts 
A' vortex) 

) 


Bulk SC 


FIG. 3: Forces on a vortex in a homolaminate. As the film is 
made thinner, the image antivortex to the right of the film has 
a stronger pull on the vortex, lowering the barrier to vortex 
penetration. 

The differing penetration depths in the layers of a het¬ 
erolaminate cause it to behave differently than a homo¬ 
laminate. Here we consider structures in which the bulk 
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has a smaller penetration depth than the film. For such 
structures, if the film is very thin, it does not provide 
much screening for the bulk, and Bi reaches the bulk’s 
Bsh before the thin film barrier disappears. As with a 
homolaminate, a very thick film behaves like a bulk, and 
reaches that material’s bulk Bsh while Bi is still relatively 
small. However, between these two extremes, there is a 
situation in which the film provides some screening, so 
that Bi is large but still smaller than Bq. In this case, a 
benefit can be realized - the small penetration depth of 
the material in the bulk causes Bi to be larger than with 
the exponential decay expected for a thick film [Eq. ([^]. 
This in turn reduces the magnitude of the negative gra¬ 
dient in Bm, bolstering the barrier to flux penetration 
(Eqn.|^. This increase in the barrier is depicted in Fig.|^ 
The dark curves show Bm, By, and Q for a NbaSn thin 
film/insulator/Nb bulk SIS structure with 10 nm thick 
insulator and d/Xf = 0.64 (the peak of the cyan curve 
in Fig. [^. The light curves show calculations for a bulk 
NbaSn film (for this case, the dark shaded region repre¬ 
senting the insulator does not apply). In this example, 
Bq = 300 mT. The Gibbs free energy of the SIS’ struc¬ 
ture is still sharply peaked, showing a relatively robust 
energy barrier, but that of the bulk film is almost flat, 
showing that flux penetration is likely to occur at slightly 
higher fields. 



FIG. 4: Comparison of a SIS’ film with near-optimal param¬ 
eters (dark curves) to a bulk film (light curves). The slower 
decay of Bm in the large-A thin film influences Q, bolstering 
the barrier to flux penetration. Note that the dark shaded 
region representing the insulating region of the SIS’ can be 
ignored for the bulk film. 


The impact of this is a modest increase in for 

these structures compared to the bulk value of the film 
material. However, the range of film thicknesses over 
which the increase is appreciable (> few %) is relatively 
small, and the gain decreases as the thickness of the in¬ 
sulating layer increases. 


IV. THIN FILMS IN THE GINZBURG-LANDAU 
APPROACH 


Calculating Bgh using London theory, as done in the 
previous section, fails to take into account 2D instabil¬ 
ities in the order parameter, therefore overestimating 
Bsh in many circumstances. The problem of calculat¬ 
ing Bsh for bulk samples while taking into account 2D 
instabilities has a long history (see, e.g., Ref. IT^ and 
has mostly been tackled in the Ginzburg-Landau (GL) 
framework. Only r ecent ly calculations beyond GL the¬ 
ory were performed p^ * ^^ * they showed that while the GL 
results cannot be trusted quantitatively at low tempera¬ 
tures, they give a qualitatively correct estimate. There¬ 
fore, for simplicity we restrict ourselves to GL theory 
even in the low-temperature regime where its quantita¬ 
tive predictions are not exact. 

The approach we use to find Bsh is described in detail 
in Ref. 1191 we first extremize the GL free energy, a func¬ 
tional of the spatially dependent order parameter ^ and 
superfluid velocity Vg, and then study its stability against 
small perturbation of these functions. In the present case, 
the GL free energy is the sum of the contributions for the 
bulk and the film. The boundary conditions are the usual 
ones for the order parameter (vanishing of its derivatives 
at all surfaces); the superfluid velocity vanishes deep into 
the bulk, and its derivative at the external film surface is 
proportional to the applied magnetic field. Similarly, the 
field between the film and the bulk is proportional to the 
derivatives of the superfluid velocities at the two surfaces. 
However, this internal field is not externally imposed, but 
must be calculated consistently with Maxwell equations; 
this gives the final condition of continuity of the vector 
potential across all surfaces. Hence for a very thin insu¬ 
lating barrier the superfluid velocities at the bulk surface 
and the internal film surface coincide, while the film su¬ 
perfluid velocity would be higher for a thicker insulator. 

For an analytical estimate of the SIS’ superheating 
field Hf/"® , we consider the simplest possible case of a 
strongly type H bulk material (GL parameter kgl ^ 1) 
shielded by an insulator of negligible thickness and a 
strongly type H thin film with <C d ^ A/. With a 
thin film, the difference between the internal field (at the 
bulk surface) and the applied field (outside the film) is 
small, and the maximum possible field is reached when 
the internal field coincides with the bulk superheating 
field Bsh,b- Indeed, within GL theory and at linear order 
in d/Xf, using the boundary conditions discussed above 
we find [see Appendix [A]: 


B 


SIS’ 


= B, 


h.b 



1 - 


^2 (a/) 


( 8 ) 


where Vs,r = is the ratio of the maximum su¬ 

perfluid velocities for bulk and film material, respectively. 
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This ratio can be written in terms of critical fields and 
penetration depths as Ug.r = BcjXf, and as a nec¬ 

essary condition for metastability it must satisfy Vg^r < 
since in the bulk material the superfluid velocity has al¬ 
ready reached its maximum possible value at the sur¬ 
face, the film material must be able to support a higher 
superfluid velocity. We stress again that for sufficiently 
thin films (below the critical thickness discussed in the 
next paragraph), as the applied field becomes larger than 
, the bulk becomes unstable, while the film is still 
(meta)stable. As qualitatively expected, Eq. Q shows 
that for better screening a thicker film should be used, 
and that as the film material penetration depth increases, 
its screening power decreases. Also, the need of small 
implies that the film material critical field should be suf¬ 
ficiently large, Bcj > Bc,bXb/Xf. Interestingly, based on 
the values reported in Table |lj this condition can be met 
if using NbaSn or MgB 2 to shield Nb. 

We note, however, that there is in principle a limit 
on how thick the film can be made: since the supercur¬ 
rent velocity at the film external surface increases with 
thickness, if the film is too thick it will become unsta¬ 
ble at a field below that predicted by Eq. Q. Within 
our approximations, we find that the critical thickness 
for the film to also become unstable at is dc = 

Xf^/^{l-Vs^r)B,j/B,^b- We see that the condition 
d < dc can severely restrict the maximum film thickness 
only in the regime Bcj Bc^b, Xf ^ A;,. For the ma¬ 
terial parameters in Table |T] our formula gives dc ~ Ay, 
but films of this thickness are beyond the approximate 
analytical treatment that leads to Eq. ([^. Therefore, to 
study the screening properties of films of intermediated 
thickness, d ^ Ay, in the next section we resort to nu¬ 
merical calculations that also account for the finite value 
of kgl- 


V. FILMS OF INTERMEDIATE THICKNESS 

For films of intermediate thickness, numerical tech¬ 
niques are needed to accurately estimate the effective 
superheating field of SIS structure in Ginzburg-Landau 
theory. Here, we follow closely the methods described in 

Ref.im 

The Ginzburg-Landau equations are solved in each 
domain separately and then boundary conditions are 
matched. In order to improve numerical stability, we im¬ 
plement the boundary conditions as follows: at the film 
surface the gradient of the order parameter is fixed to 
zero while the magnitude is allowed to vary, effectively 
defining the applied magnetic field implicitely in terms 
of the order parameter. We also allow the value of the 
order parameter and the vector potential on the film side 
of the interface to vary. On the bulk side of the interface, 
the gradient of the magnetic field is fixed to zero while 
its magnitude is allowed to vary. Infinitely deep in the 
bulk the order parameter is fixed to one and the vector 
potential vanishes. This configuration introduces three 


parameters for the boundary conditions: the magnitude 
of the order parameter on either side of the interface and 
the magnitude of the vector potential on the film side of 
the interface. These three parameters are varied until the 
gradient of the order parameter vanishes on the film side 
of the interface and both the magnetic vector potential 
and the magnetic field are continuous at the interface. 

Having found a solution, we next solve the eigenvalue 
problem associated with the stability of the solution to 
infinitesimal fluctuations of wavenumber k as in Ref. [IS 
These solutions are also found numerically using bound¬ 
ary conditions similar to those just described. The mag¬ 
nitude of the applied magnetic field and the wavenumber 
are then varied simultaneously to identify the least stable 
fluctuation and the applied magnetic field at which it just 
becomes unstable (i.e., at which the eigenvalue becomes 
zero). In this way we identify the superheating field and 
the critical wavenumber that characterizes the unstable 
fluctuations. These calculations are summarized in Fig.[^ 
in which we plot as a function of film thickness for 

various materials. Note that the dashed lines start from 
thicknesses of about 50 nm. For films with thickness less 
than this, numerical results become increasingly difficult, 
presumably due to the extremely separated length scales 
involved. Interestingly this thickness coincides with ap¬ 
proximately twice the depth ^Xf^f of the fluctuation^, 
suggesting that interactions between the fluctuations of 
both films surfaces may become relevant. Moreover, nu¬ 
merical solutions indicate that at finite k, the nature of 
the of the instability itself may change from 2D to ID as 
the thickness decreases. Although beyond the scope of 
the present work, these indications deserve further inves¬ 
tigation. 
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FIG. 5: Bsh in a SIS’ structure as a function of film thickness. 
The insulating layer is assumed to be very thin. London limit 
calculations are compared to Ginzburg-Landau analytical and 
numerical calculations. For reference, the dashed vertical line 
is at position d ~ 5^/ ~ y/Xf^f. We remind that the results 
of Secs. |IIl| (solid lines) and |IV| (dot-dashed line) are valid for 
films thick compared to 

The Ginzburg-Landau calculations show good qualita¬ 
tive agreement with the London calculations from Fig. 
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also shown in this figure. There are some quantitative 
differences, likely due to the approximations used in the 
London limit. For instance, the difference in the calcu¬ 
lated bulk Bsh of the film material, which is approached 
as the him becomes a few A/ thick, is due to hniteness 
of K. For the heterolaminate, in both cases as the him 
thickness increases, Sf/"® shows a peak near d ~ A/, 
then decreases to the superheating held of the him ma¬ 
terial as the him becomes very thick. The thickness at 
which the peak occurs is somewhat smaller for the Lon¬ 
don limit, but the two plots are otherwise very similar in 
shape. 


VI. CONCLUSIONS 

In this study, we analyzed the magnetic shielding prop¬ 
erties of superconductors at high helds and high fre¬ 
quencies. To prevent strong vortex dissipation due to 
drag, the analysis was restricted to a regime where hux 
penetration is not allowed. The London-limit numeri¬ 
cal results were verihed against analytical and numeri¬ 
cal Ginzburg-Landau calculations. We showed that the 
SIS’ structure can produce a modest enhancement of the 
maximum screening field compared to a single supercon¬ 
ducting slab for certain materials and film thicknesses, 
see the maxima in Fig. 
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Appendix A: Derivation of Eq. [8| 


take the insulator thickness to be negligible, (5 = 0, as this 
gives the highest possible metastable field. As discussed 
in Secs. |T^ and |TVj the instability happens at the bulk 
surface; this fixes the values of the superfluid velocity at 
the interior surface of the film to be: 




^c,b^b 


(A2) 


Clearly, a necessary condition for the metastability of 
the him is Vs^r < 1- In fact, since the superhuid velocity 
at the outer surface is larger, we will further need to 
check that < 7 g(—d/2A/) < 1/3. In addition to the above 
boundary condition, we also need the held between him 
and bulk to coincide with the bulk superheating held: 


/ f ^ \ _ ^sh,b 


(A3) 


The task is now to hnd the external held at which these 
two boundary conditions are satished, 



(A4) 


To solve Eq. AT thanks to the assumption d <C A/, we 
can proceed by a Taylor expansion of the function qo{x) 
near x = Q: 


Qo{x) = Qc + box + bi— + b2— -I- ... (A-5) 

Substituting the expansion into Eq. |Al| and matching the 
terms on the two side of the equality we hnd 


bi = gc (l - ql) , b2 = (l - 3g/) , (A6) 


showing that only two parameters of the expansion, q,. 
and bo, are left undetermined and thus can be hxed by the 
boundary conditions. Moreover, Eq. |A4| can be written 
in the form 


In the case of large GL parameter kql ^ 1, the cal¬ 
culation of the metastability held, Eq. is greatly sim- 
plihed: for kql ^ oo, the spatial prohle of the order pa¬ 
rameter is fully determined by that of the superhuid ve¬ 
locity, and the differential equation for the latter is local, 
albeit non-linear.^i^ Indicating with qo the dimensionless 
superhuid velocity, for the geometry we are considering 
it obeys the equation: 


tdSIS' _ p 

^sh — ^sh,b 


-V2B,jq,(l-ql)±+o(^^^ , (A7) 


and hence to calculate Rf/"® to second order in d/A/ we 
only need to know q^ to hrst order. We can therefore use 
the boundary condition (A31 at lowest order to obtain 


bo = Bsh,b/V‘^Bcj and the boundary condition (A2) at 
hrst order to hnd 


q'd = 90-90 (Al) 

and the metastability condition takes the simple form 
9 q < 1/3.1^ The dimensionful velocity is proportional 
to 9 o multiplied by the critical held and the penetration 
depth, Vs oc B^X, and we will not need the proportional¬ 
ity constant in what follows. 

For simplicity, in this Appendix we use a coordinate 
system in which x-axis perpendicular to the him has its 
origin in the middle of the him and measure lengths in 
units of the him material penetration depth A/. We also 


9c = - 



V2B^j ^ ■ 


(A8) 


Substituting this expression into Eq. |A7| and keeping only 
terms up to second order we hnd 


B 


SIS’ 


= B, 


h.b 


1+^(1- <c) 







d 


(A9) 
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To put this equation in the form given in Eq. we qo{—d/2Xf) > —1/-\/3, we obtain the critical thickness 

perheating and critical fields. Finally, by considering 
at linear order in d/Xf the metastability requirement 


use the relationship!!^ Bsh,b = v5Bc,b/^ between su- dc reported at the end of Sec. IV 
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With flux penetration field we mean the maximum value 
of the magnetic field just outside the film such that no flux 
enters into the bulk. 

Niobium has fairly low kappa, so the London limit does 
not apply. In these calculations, the limit used for bulk 
niobium is when the surface reaches Bsh from Table |T] For 
all other materials, the slope of the Gibbs free energy is 
used as the criterion for flux penetration as described. 






